Orthogonality Catastrophe for Vortices in rf-Wave Superconductors 
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The dynamics of magnetic vortices in the mixed state of d-wave superconductors is affected 
by interaction with quasiparticles near the gap nodes. We study this effect by computing the 
overlap of the ground state wave functions of the nodal quasiparticles in a two-dimensional d-wave 
superconductor as the ground state changes in response to the motion of the vortex. We find that 
the overlap is strongly suppressed. This orthogonality catastrophe is specific to a half-quantum 
vortex and disappears in the case of a doubly-quantized vortex. This implies strong inhibition of 
the motion of half-quantum vortices. The results allow us to develop a simple description of the 
quantum vortex creep at low temperatures. 

PACS numbers: 74.25.Qt 



Magnetic vortices in Cuprate superconductors 0, ||] 
behave rather differently compared to their counterparts 
in conventional materials. One of the characteristic fea- 
tures of the superconducting state of the Cuprates is the 
existence of the low-energy quasiparticles near the nodes 
of the d-wave gap. One is naturally led to ask: what is 
the influence of the interaction between the vortices and 
nodal quasiparticles on vortex dynamics? 

Before answering this question, let us briefly describe 
our current knowledge. Vortex creep at low tempera- 
tures is believed to proceed by quantum tunneling. It 
has been argued that two forces govern vortex dynamics 
during tunneling: the Magnus force j^] and the friction 
force. Their ratio determines the "cleanness" of the ma- 
terial. The creep in the dirty, Q super-clean, and 
intermediate regimes has been studied theoretically. 

Quasiparticles play prominent role in the vortex dy- 
namics. In s-wave superconductors, finite gap in the bulk 
precludes the existence of low-energy extended states, 
however quasiparticle states localized in the cores of vor- 
tices, where the gap vanishes, exist. [Q Transitions be- 
tween these states give rise to the friction force and 
the renormalization of the Magnus force. [|| The sit- 
uation with quasiparticle states localized inside vortex 
cores in d-wave superconductors is far less clear. ^ On 
the other hand, extended states exist and are expected 
to strongly interact with vortices. Indeed, they have 
been argued |ll| to give a quantized value of ther- 
mal Hall conductance in the mixed state at low tem- 
peratures. This prediction is hard to verify since current 
experiments are performed at relatively high temper- 
atures. Quantum vortex creep, which is also expected 
to be sensitive to the interaction with quasiparticles, oc- 
curs at much lower temperatures and thus offers experi- 
mental information about the quasiparticle spectrum at 
a much lower energy scale. 

To study the contribution of nodal quasiparticles to 
vortex action we calculate the overlap of the ground state 
wave functions of the nodal fermions as the ground state 
changes in response to the motion of the vortex. We 



find that this overlap is extremely small in the case of 
a half-quantum vortex. This implies that the motion 
of half-quantum vortices is strongly suppressed due to 
the interaction with nodal quasiparticles. Indeed, we 
argue that this orthogonality catastrophe provides the 
dominant contribution to the vortex action at tempera- 
tures relevant to quantum vortex creep. Based on this, 
we present a simple phenomenological description of the 
available experimental data. 

Let us consider a vortex in a 2D d-wave superconduc- 
tor. We use the Bogoliubov-de Gennes (BdG) equation 
to treat the interaction of quasiparticles with the vor- 
tex. ||l^, Q . We first take advantage of the conservation 
of spin in the BCS theory by changing from the original 
electronic operators c to the new fermions: V'lfe = Cfct, 
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The new operators are combined to form 



a two-component spinor: ^ = (■0i,V'2)- The density of 
the '0-particles t/'^V' is the z-component of the spin den- 
sity of the original electrons, and visa versa. Next, we 
concentrate on the variation of the phase of the order 
parameter 4>{r) around the vortex and neglect the vari- 
ation of the amplitude. We then eliminate the phase 
winding by making a single-valued gauge transformation: 
-01 exp[— i(/)(r)]i/)i, V2 V'2- The BdG equation takes 
the form Htp = eip with the following Hamiltonian (we 
use units with h = c = 1): 



H = E{p + a + js(Jz)<7z + A(p + d)cr^. 



(1) 



Here, az , (7^ are Pauli matrices, p — (px ,Py) is the 2D mo- 
mentum, E{p) is the bare electron energy relative to the 
Fermi energy and A(p) is the gap function. The presence 
of the vortex produces two contributions to the Hamilto- 
nian: 1) There appears a "Doppler shift". The superflow 
js = 90/2 — eA shifts the energies of the nodal quasipar- 
ticles similar to the classical Doppler effect. 2) The other 
contribution is of topological nature. It requires that 
a parallel transport of a quasiparticle around the vor- 
tex should lead to the change of its phase by tti/, where 
v is the vortex quantum number. This is enforced by 
the Chern-Simons gauge field a — 90/2 which satisfies 
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d X a = 7rz^M(r) (assuming the vortex is at the origin). 

First, we would Hke to calculate the overlap between 
the ground state wave functions of the superconductor 
with and without the half-quantum vortex. The two 
terms in the Hamiltonian that were discussed above both 
perturb the system shifting the single quasiparticle eigen- 
states. An important difference between them is their 
spatial range. The superflow js(r) is screened at dis- 
tances larger than the penetration depth. On the con- 
trary, the field a{r) decays as 1/r at all distances. Let 
us consider the Doppler shift term first. The solution of 
the orthogonality catastrophe problem in the case of 
a short-range perturbing potential is given by: 



l(/K>l 

A 



(2) 
(3) 



Here, \i) and |/) are the ground states of the system 
before and after the perturbing potential is switched on, 
N is the total number of fcrmions in the system, and 
S{e) — (l/2z) lnS'(£) is the phase shift matrix defined as 
the logarithm of the scattering matrix. The trace in the 
Eq. (||) is over all states at the Fermi energy. In our case 
the Fermi energy is zero. The density of states of nodal 
quasiparticles vanishes linearly with energy: there are no 
states at zero energy. This property of the phase space 
implies that the orthogonality exponent A is zero for any 
scattering potential of finite range. This means that the 
Doppler shift term causes possibly quite large but finite 
reduction of the overlap. Absence of infrared divergence 
in this case has been found previously in a different 
context. 

Based on what was said above, we set in the following 
js = in Eq. (^. [l4| The remaining perturbation of 
the system is the long-range vector Chern-Simons field 
a(r). Because the field decays slowly at infinity, it is 
not clear whether the above result for the orthogonality 
catastrophe for short-range potentials can be generalized 
to this case. It is also not clear whether at all one can 
use plane waves as asymptotic states to calculate the S 
matrix: the asymptotic completeness has been proven 
only for short-range potentials in the Dirac equation. 

For these reasons we shall adopt another approach. 
The field a(r^ can be gauged away by a gauge transforma- 
tion: ip t/' 6xp(w0/2), where v is the vortex quantum 
number and 9 is the polar angle. The BdG Hamiltonian 
Eq. (|l|) reduces to a free one: 



(4) 



For a doubly-quantized vortex, v — 2, this implies no or- 
thogonality catastrophe between the ground states with 
and without such a vortex. On the other hand, for a 
singly-quantized vortex this transformation makes ip non- 
single valued. To take this into account, we introduce a 
cut C in the {x,y) plane as shown in Fig. 1. The vor- 
tex is assumed to be at point A at the origin. Since the 
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FIG. 1: The cut C. 

cut should end on another half-quantum vortex, we in- 
troduce one at point B for definiteness. We then impose 
a TT-phase shift across the cut. 

The distance L between the two vortices can be in- 
terpreted in various ways: First, L may be thought of 
as a long-distance cut-off. In another setting, one may 
consider an isolated vortex being moved from point A to 
point B. This can be done by removing a vortex at point 
A (inserting a vortex with negative vorticity) and creat- 
ing one at point B. In this situation L is the distance by 
which the vortex has been moved. 

Having made the gauge transformation we can refor- 
mulate our original problem as that of finding the over- 
lap of the ground state wave functions of the system be- 
fore and after a 7r-phase shift across the cut shown in 
Fig. 1 was imposed. To this end, consider the dynami- 
cal overlap function which is defined as follows: Let \i) 
be the ground state wave function of the unperturbed 
system. We shall work in imaginary time to simplify 
calculations and avoid problems with time-ordering. At 
time r = the phase shift tt across the cut C is im- 
posed. Let |t) — C/(r)|i), where U{t) is the imaginary 
time evolution operator. Since in the limit r — > cxd only 
the contribution from the new ground state to U (r) sur- 
vives, we have: limT-^oo('''|i) = (/I*): where |/) is the 
new ground state. Since (r|i) = (i|[/t(T)|i), the problem 
reduces to that of finding the expectation value of the 
evolution operator in the old ground state. We have: 



U{t) = exp{i7r(5(T)} , 



Q{r) 



J{T')dT'. 



(5) 
(6) 



Here, J(t) is the total current through the cut C and 
Q{t) is the net charge that has fiowed through the cut 
between times and r (here and below, unless mentioned 
otherwise, the charges and currents are those of the "tp- 
fermions). The meaning of the above expressions is quite 
transparent: Each fcrmion contributes phase tt to the 
path integral every time it crosses the cut. The opera- 
tor U{t) can thus be thought as a generator of counting 
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statistics, pq ] 

Since the unperturbed Hamiltonian is quadratic in the 
fermionic currents we can apply Gaussian averaging to 
obtain: 



{T\i) = cxp{-v4}, 



(7) 
(8) 



To proceed, we need to find the current-current corre- 
lation function. We linearize the Hamiltonian Eq. (^) 
around one of the nodes (e.g. one at {kp,{))): E{p) = 
vfPx, ^(p) = VAPy The anisotropy parameter a = 
vf/va ~ -Bf/A is quite large (a w 14 for YBCO). We 
then make a = 1 by rescaling the distances: Px Px/vp, 
Py —^Py/vA- This allows us to exploit the full "relativis- 
tic" invariance of the problem. Indeed, the imaginary 
time action (in the second-quantized notation) becomes: 



(9) 



Here, = {w,Px,Py) is the 2+1 momentum, ip = 
—iip'^cTy. Pauli matrices cr^ form a Clifford algebra: 
{a^,a^} — 25 ^y. The "rclativistic" invariance is now 
explicit since both cr^ and Pf^ transform as 2+1 vectors 
under rotations of the Euclidean space-time and the ac- 
tion S is therefore a scalar. The conserved current is: 



3^t = i^cr^tp. 



(10) 



The current-current correlator is given (in the free the- 
ory) by the so-called bubble or vacuum polarization dia- 
gram: 



(2^ 



Tr [a^G{q)a,G{q + k)]. (11) 



Here, G{q) = (T^f/p/q^ is the Green function of the Dirac 
equation. The correlator, as written, is given by an ultra- 
violet divergent expression. This divergence can be elim- 
inated by an appropriate high-energy cut-off. Introduc- 
tion of such a cut-off, however, will break the gauge in- 
variance. A more physical way is to explicitly impose the 
current conservation condition: k^{jfj_j^)k = 0. This way 
one obtains the well-known result (see for a related 
calculation) : 



The inverse Fourier transform of this correlator is: 

1 



{J^.{r)J.m = ~{^'5,,~^^^,) 
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(12) 



(13) 



where — (t',x,?/). The fluctuation of the charge 
{Q'^{t)) that flowed through the cut C of length L dur- 
ing time r is then given by the integral over the segment 



r2 = {r : < r' < r, a; = 0, < y < L} of the hyperplane 
{T,y). Explicitly, one has: 



dri / dr2 {jxiri)jx{r2)), 



(14) 



Substituting Eq. ( |13D into Eq. ( |14| ) and using Gauss the- 
orem twice, one obtains: 



1 
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dn Jdn 



dS{n)dS{r2) 
(ri - r2)2 



(15) 



Here, dS = Adl is the "surface" element of the boundary 
90, while dl is the length element of the corresponding 
perimeter. The integral in Eq. (|l5| ) diverges as ri r2- 
This means that the exponent in the final result acquires 
an explicit cut-off dependence. A suitable way to intro- 
duce a short time cut-off is to replace: 



1 



1 



(16) 



The cut-off a is, roughly, the inverse of the energy at 
which the linear Dirac dispersion breaks down. It would 
be justified to assume a ~ 1/A. The integral in Eq. ( |T5| ) 
can now be easily evaluated. Plugging the result in 
Eq. (^) in the case where L, r ;3> A^^ one obtains: 



Ar^A{L + T) 



(17) 



This can be called a perimeter law for vortex action: 
the contribution of nodal fermions to the imaginary time 
action of a half-quantum vortex is proportional to the 
length of the vortex trajectory in the Euclidean space- 
time. The proportionality of the action to the perimeter 
of n rather than its surface area is essentially due to 
the conservation of the number of V'-particles. In terms 
of the original electrons, it is the conservation of the z- 
component of spin polarization. 

Let us now discuss the implication of these results for 
the vortex dynamics in d-wave superconductors. First 
of all, Eq. (17) shows that the orthogonality catastrophe 
for a half-quantum vortex is very strong. Even when the 
length L is finite (e.g. when the vortex moves over a 
small distance) the dynamical overlap of the two states 
vanishes in the limit r — » oo. Certainly, finite temper- 
ature provides a cut-off: nh ^ l/T. This leads to the 
conclusion that the contribution of the nodal fermions to 
the imaginary time action at low temperatures is domi- 
nated by the second term in Eq. (^: A ^ A/T. This, 
in turn, means that the dynamical magnetization relax- 
ation rate Q{T) = l/S, ^ where S is the total tunneling 
action, vanishes linearly with temperature in the limit 
T ^ 0. Indeed, the diverging contribution A to the ac- 
tion S becomes the dominant one in the same limit. The 
fact that experiments ||^ show a finite <5(0) leads us to 
conclude that there is some other low-energy cut-off be- 
sides temperature. We speculate that it may appear as 
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FIG. 2: Qualitative behavior of the dynamical magnetization 
relaxation rate Q as a function of temperature T. 



scenario is similar to the one used in recent studies of frac- 
tionalized phases and deserves further investigation. 

In summary, we have found that the interaction be- 
tween magnetic vortices and nodal quasiparticles in d- 
wave superconductors leads to a strong orthogonality 
catastrophe in response to the motion of half-quantum 
vortices. This effect strongly suppresses their mobil- 
ity and is argued to give the dominant contribution to 
the vortex action at low temperatures. Based on these 
results, we have developed a simple description of the 
quantum creep of vortices in d-wave superconductors. 
We would like to thank A. Paramekanti, A. Ludwig and 
M. P. A. Fisher for discussions and criticism. 



a consequence of a mini-gap in the Dirac spectrum of 
nodal quasiparticles. Such a mini-gap is expected in fi- 
nite magnetic field as a result of curvature terms in the 
dispersion around the nodal point. Its magnitude has 
been estimated ||T^ to be ^ = kH, k « O.SiiTTesla"^, 
where H is the external magnetic field. Taking the mini- 
gap into account, we expect the following behavior of the 
magnetization relaxation rate at low temperatures: 



QiT) 



r T/A, 



T>5 
T<5 



(18) 



This behavior is shown in Fig. 2. Qualitatively, the 
following picture emerges: There is a small energy scale 
6 set by the size of the mini-gap in the Dirac spectrum 
of nodal fermions. This energy scale defines a cross-over 
temperature between two regimes: a) For T ^ S the 
function Q{T) is linear in T with the coefficient ^ 1/A 
which is little dependent on the details; b) For T ^ S 
the function Q{T) saturates at a value ~ 6 /A. Experi- 
mentally, [|| in the regime T ^ 6 typical slopes of Q{T) 
are lO^^K^^ which is consistent with the typical value 
A ~ lO^K for the Cuprates, although the linearity in T 
is difficult to establish. Typical values of the cross-over 
temperature are 6 ~ IK (i.e. for magnetic fields ~ IT, in 
agreement with the theoretical estimate 1 10 ) . The argu- 
ments above then predict Q{Q) ~ 10~^, also in agreement 
with experiments. Note that in this simple picture we 
have omitted disorder effects on the spectrum of nodal 
quasiparticles, inter-vortex interactions etc. In particu- 
lar, the Zeeman splitting is expected to change the 
simple S = kH behavior. 

As an end note, we would like to point out one inter- 
esting consequence of the obtained results. As we have 
seen, the orthogonality catastrophe strongly inhibits the 
motion of half-quantum vortices, whereas the effect dis- 
appears for doubly-quantized vortices. This suggests a 
possible tendency of vortices to pair which would increase 
their mobility and, effectively, decrease their energy. This 
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